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Abstract. We introduce an extension to Merton's famous continuous time model of optimal 
consumption and investment, in the spirit of previous works by Pliska and Ye, to allow for 
a wage earner to have a random lifetime and to use a portion of the income to purchase life 
insurance in order to provide for his estate, while investing his savings in a financial market 
comprised of one risk-free security and an arbitrary number of risky securities driven by multi- 
dimensional Brownian motion. We then provide a detailed analysis of the optimal consumption, 
investment, and insurance purchase strategies for the wage earner whose goal is to maximize 
the expected utility obtained from his family consumption, from the size of the estate in the 
event of premature death, and from the size of the estate at the time of retirement. We use 
dynamic programming methods to obtain explicit solutions for the case of discounted constant 
relative risk aversion utility functions and describe new analytical results which are presented 
together with the corresponding economic interpretations. 



1. Introduction 

We consider the problem faced by a wage earner having to make decisions continuously about 
three strategies: consumption, investment and life insurance purchase during a given interval 
of time [0, min{T, r}], where T is a fixed point in the future that we will consider to be the 
retirement time of the wage earner and r is a random variable representing the wage earner's 
time of death. We assume that the wage earner receives his income at a continuous rate i{t) 
and that this income is terminated when the wage earner dies or retires, whichever happens 
first. One of our key assumptions is that the wage earner's lifetime r is a random variable and, 
therefore, the wage earner needs to buy life insurance to protect his family for the eventuality 
of premature death. The life insurance depends on a insurance premium payment rate p{t) such 
that if the insured pays p{t) ■ 6t and dies during the ensuing short time interval of length 6t 
then the insurance company will pay p{t)/rj[t) dollars to the insured's estate, where 'q{t) is an 
amount set in advance by the insurance company. Hence this is like term insurance with an 
infinitesimal term. We also assume that the wage earner wants to maximize the satisfaction 
obtained from a consumption process with rate c{t). In addition to consumption and purchase 
of a life insurance policy, we assume that the wage earner invests the full amount of his savings 
in a financial market consisting of one risk-free security and a fixed number > 1 of risky 
securities with diffusive terms driven by M-dimensional Brownian motion. 

The wage earner is then faced with the problem of finding strategies that maximize the utility 
of (i) his family consumption for all t < minjT, r}; (ii) his wealth at retirement date T if 
he lives that long; and (iii) the value of his estate in the event of premature death. Various 
quantitative models have been proposed to model and analyze this kind of problem, at least 
problems having at least one of these three objectives. This literature is perhaps highlighted by 
Yarri [9] who considered the problem of optimal financial planning decisions for an individual 
with an uncertain lifetime as well as by Merton [H [5] who emphasized optimal consumption and 
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investment decisions but did not consider life insurance. These two approaches were combined 
by Richard, whose impressive paper [8] uses sophisticated methods at an early date for the 
analysis of a life-cycle life insurance and consumption-investment problem in a continuous time 
model. Later, Pliska and Ye [Ul [7] introduced a continuous-time model that combined the 
more realistic features of all those in the existing literature and extended the model proposed 
previously by Richard, the main difference being the choice of the boundary condition, leading to 
somewhat different economic interpretations of the underlying problem. More precisely, while 
Richard assumed that the lifetime of the wage earner is limited by some fixed number, the 
model introduced by Pliska and Ye had the feature that the duration of life is a random variable 
which takes values in the interval ]0,oo[ and is independent of the stochastic process defining 
the underlying financial market. Moreover, Pliska and Ye made the following refinements to the 
theory: (i) the planning horizon T is now seen as the moment when the wage earner retires, 
contrary to Richard's interpretation as a finite upper bound on the lifetime; and (ii) the utility of 
the wage earner's wealth at the planning horizon T is taken into account as well as the utility of 
lifetime consumption and the utility of the bequest in the event of premature death. Blanchet- 
Scalliet et al. paper [1] deals with optimal portfolio selection with an uncertainty exit time for a 
suitable extension of the familiar optimal consumption investment problem of Merton, without 
considering any kind of life insurance purchase. 

Whereas Pliska and Ye's financial market involved only one security that was risky, in the 
present paper we study the extension where there is an arbitrary (but finite) number of risky 
securities. The existence of these extra risky securities gives greater freedom for the wage earner 
to manage the interaction between his life insurance policies and the portfolio containing his 
savings invested in the financial market. Some examples of these interactions are described 
below. 

Following Pliska and Ye, we use the model of uncertain life found in reliability theory, com- 
monly used for industrial life-testing and actuarial science, to model the uncertain time of 
death for the wage earner. This enables us to replace Richard's assumption that lifetimes are 
bounded with the assumption that lifetimes take values in the interval ]0, oo[. We then set up 
the wage earner's objective functional depending on a random horizon min{T, r} and trans- 
form it to an equivalent problem having a fixed planning horizon, that is, the wage earner who 
faces unpredictable death acts as if he will live until some time T, but with a subjective rate 
of time preferences equal to his "force of mortality" for his consumption and terminal wealth. 
This transformation to a fixed planning horizon enables us to state the dynamic programming 
principle and derive an associated Hamilton-Jacobi-Bellman (HJB) equation. We use the HJB 
equation to derive the optimal feedback control, that is, the optimal insurance, portfolio and 
consumption strategies. Furthermore, we obtain explicit solutions for the family of discounted 
Constant Relative Risk Aversion (CRRA) utilities and examine the economic implications of 
such solutions. 

In the case of discounted CRRA utilities our results generalize those obtained previously by 
Pliska and Ye. For instance, we obtain: (i) an economically reasonable description for the optimal 
expenditure for insurance as a decreasing function of the wage earner's overall wealth and a 
unimodal function of age, reaching a maximum at an intermediate age; (ii) a more controversial 
conclusion that possibly an optimal solution calls for the wage earner to sell a life insurance 
policy on his own life late in his career. Nonetheless, the extra risky securities in our model 
introduce novel features to the wage earner's portfolio and insurance management interaction 
such as: (i) a young wage earner with small wealth has an optimal portfolio with larger values 
of volatility and higher expected returns, with the possibility of having short positions in lower 
yielding securities; and (ii) a wage earner who can buy life insurance policies will choose a more 
conservative portfolio than a wage earner who is without the opportunity to buy life insurance, 
the distinction being clearer for young wage earners with low wealth. 
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This paper is organized as follows. In section[2]we describe the problem we propose to address. 
Namely, we introduce the underlying financial and insurance markets, as well as the problem 
formulation from the point of view of optimal control. In section [3] we see how to use the 
dynamic programming principle to reduce the optimal control of section [2] to one with a fixed 
planning horizon and then derive an associated HJB equation. We devote section H] to the case 
of discounted CRRA utilities. We conclude in section [SI 

2. Problem formulation 

Throughout this section, we define the setting in which the wage earner has to make his 
decisions regarding consumption, investment and life insurance purchase. Namely, we introduce 
the specifications regarding the financial and insurance markets available to the wage earner. 
We start by the financial market description, followed by the insurance market and conclude 
with the definition of a wealth process for the wage earner. 

2.1. The financial market model. We consider a financial market consisting of one risk- 
free asset and several risky-assets. Their respective prices (5'o(t))o<t<T and iSn{t))o<t<T for 
n = 1, N evolve according to the equations: 

dSoit) = r{t)So{t)dt , 5o(0) = so , 

M 

dSn{t) = lln{t)Sn{t)dt + Sn{t) a„^(t)dW^(t) , 5„(0) = S„ > , 

m=l 

where W{t) = {Wi{t), . . . , VFA/(t))^ is a standard M-dimensional Brownian motion on a prob- 
ability space {Q,,T,P), r{t) is the riskless interest rate, = {fii{t), . . . , fj,]\f{t)) G is the 
vector of the risky-assets appreciation rates and a{t) = {(Tnm{t))i<n<N,i<m<M is the matrix of 
risky-assets volatilities. 

We assume that the coefficients r{t), ^{t) and a{t) are deterministic continuous functions on 
the interval [0, T\. We also assume that the interest rate r{t) is positive for all t G [0, T] and the 
matrix a{t) is such that aa'^ is nonsingular for Lebesgue almost all t G [0, T] and satisfies the 
following integrability condition 

N M „r 
n=lm=l ''^ 

Furthermore, we suppose that there exists an (J"j)o<t<T-pi'ogressively measurable process 7r(t) G 
M*^, called the market price of risk, such that for Lebesgue-almost-every t G [0,T] the risk 
premium 

a{t) = {fii{t)-r{t),...,fiN{t)-r{t)) eR^ (1) 
is related to 7r(t) by the equation 

a{t) = a{t)TT{t) a.s. 
and is such that the following two conditions hold 

T 

llvrft)!!^ < oo a.s. 

1 „ M 

1 . 







E 



exp(^-J^ Tr{s)dW{s) ||7r(s)fds^ 



The existence of such a process 7r(t) ensures the absence of arbitrage opportunities in the financial 
market defined above. Note also that the conditions on the matrix a above do not imply market 
completeness. See ^ for further details on market viability and completeness. 
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Moreover, throughout the paper we will assume that {fl,T,P) is a filtered probability space 
and that its filtration F = {J-t,t € [0, T]} is the P-augmentation of the filtration generated by 
the Brownian motion W{t), a{W{s),s < t} for t > 0. Each sub-a-algebra J-'t represents the 
information available to any given agent observing the financial market until time t. 

2.2. The life insurance market model. We assume that the wage earner is alive at time 
t = and that his lifetime is a non-negative random variable r defined on the probability space 
{Q,J-,P). Furthermore, we assume that the random variable r is independent of the filtration 
F and has a distribution function F : [0, oo) — >■ [0, 1] with density / : [0, oo) M+ so that 



F{t)= ff{s)ds 
Jo 



Wc define the survivor function F : [0, oo) [0, 1] as the probability for the wage earner to 
survive at least until time t, i.e. 

F{t) = P{T>t) = 1- F{t) . 

Wc shall make use of the hazard function, the conditional, instantaneous death rate for the wage 
earner surviving to time t , that is 

Throughout the paper, we will suppose that the hazard function A : [0, oo) — ^ M"*" is a continuous 
and deterministic function such that 



poo 

/ X{t) di = oo 
Jo 



These two concepts introduced above are standard in the context of reliability theory and actu- 
arial science. In our case, such concepts enable us to consider an optimal control problem with 
a stochastic planning horizon and restate it as one with a fixed horizon. 

Due the uncertainty concerning his lifetime, the wage earner buys life insurance to protect 
his family for the cvcntTiality of premature death. The life insurance is available contimioTisly 
and the wage earner buys it by paying a insurance premium payment rate p{t) to the insurance 
company. The insurance contract is like term insurance, with an infinitesimally small term. 
If the wage earner dies at time r < T while buying insurance at the rate p{t), the insurance 
company pays an amount p{t) /rj^r) to his estate, where rj : [0, T] — t- M"*" is a continuous and 
deterministic function which we call the insurance premium-payout ratio and is regarded as fixed 
by the insurance company. The contract ends when the wage earner dies or achieves retirement 
age, whichever happens first. Therefore, the wage earner's total legacy to his estate in the event 
of a premature death at time r < T is given by 

where X(t) denotes the wage earner's savings at time t. 

2.3. The wealth process. We assume that the wage earner receives an income i{t) at a con- 
tinuous rate during the period [0, minjT, r}], i.e., the income will be terminated either by his 
death or his retirement, whichever happens first. Furthermore, we assume that i : [0, T] R+ 
is a deterministic Borel-measurable function satisfying the integrability condition 



/ i{t) dt < oo . 
Jo 
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The consumption process (c(t))o<t<T is a (J-i)o<t<T-pi'ogressively measurable nonnegative 
process satisfying the fohowing integrabihty condition for the investment horizon T > 

rT 



/ c{t) dt < OO a.s. . 
Jo 



We assume also that the insurance premium payment rate ip(t))o<t<T is a (J"j)o<j<r-predictable 
process, i.e., p{t) is measurable with respect to the smallest cr-algebra on x il. such that all 
left-continuous and adapted processes are measurable. In a intuitive manner, a predictable 
process can be described as such that its values are "known" just in advance of time. 

For each n = 0,1,..., and t £ [0,7"], let 9nit) denote the fraction of the wage earner's 
wealth allocated to the asset Sn at time t. The portfolio process is then given by @{t) = 
{eo{t),ei{t),--- ,0^(0) where 

N 

Y^en{t) = i, 0<t<T. (3) 

n=0 

We assume that the portfolio process is (J"j)o<t<T-pi'ogressively measurable and that, for the 
fixed investment horizon T > 0, we have that 

T 

\\Q{t)f dt < oo a.s. , 





where ||-|| denotes the Euclidean norm in 

The wealth process X{t), t E [0, min{T, r}], is then defined by 

X{t) =x+ f [^{s) - c{s) - p{s)] ds + Y, f ^"gl^f^ dSnis) , (4) 

where x is the wage earner's initial wealth. This last equation can be rewritten in the differential 
form 

dX{t) = (^i{t) - c{t) - p{t) + {9o{t)r{t) + 0n{t)fin{t)^ X{t)j dt 

N M 

+ J2 Gn{t)X{t) J2 ^nm{t)dW^{t) , (5) 
n=l m=l 

where < t < min{r, T}. 

Using relation ^ we can always write 9Q{t) in terms of . . . , 0Ar(t), so from now on we 
will define the portfolio process in terms of the reduced portfolio process 9{t) e given by 

9{t) = {ei{t),92{t),--- ,0jv(t))€M^ . 

2.4. The optimal control problem. The wage earner is faced with the problem of finding 
strategies that maximize the expected utility obtained from: 

(a) his family consumption for all t < minjT, r}; 

(b) his wealth at retirement date T if he lives that long; 

(c) the value of his estate in the event of premature death. 

This problem can be formulated by means of optimal control theory: the wage earner's goal is 
to maximize some cost functional subject to (i) the (stochastic) dynamics of the state variable, 
i.e., the dynamics of the wealth process X{t) given by ([U; (ii) constraints on the control variables, 
i.e., the consumption process c(t), the premium insurance rate p{t) and the portfolio process 
0{t); and (iii) boundary conditions on the state variables. 
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Let US denote by A{x) the set of all admissible decision strategies, i.e., all admissible choices 
for the control variables u = {c,p,9) G R-^"*"^. The dependence of A{x) on x denotes the 
restriction imposed on the wealth process by the boundary condition ^(0) = x. In particular, 
A{x) must be such that for each G A{x) the corresponding wealth process satisfies X{t) > 
for all t < min{r, r}. 

The wage earner's problem can then be restated as follows: find a strategy v = {c,p, 6) G A{x) 
which maximizes the expected utility 



V{x) = sup Eq^x 



/ U{c{s), s) ds + B{Z{t),t)I{,^t} + W{X{T))I{,^r} , (6) 
Jo 



where T At = min{T, r}, 1^ denotes the indicator function of event A, U{c, •) is the utility 
function describing the wage earner's family preferences regarding consumption in the time 
interval [0, minjT, r}], B{Z, •) is the utility function for the size of the wage earners's legacy in 
case T < T, and is the utility function for the terminal wealth at time t = T in the case 

r > T. 

We suppose that U and B are strictly concave on their first variable and that W is strictly 
concave on its sole variable. In section U] we specialize our analysis to the case where the wage 
earner's preferences are described by discounted CRRA utility functions. 

3. Stochastic optimal control 

In this section we use the techniques in [6, '2] to restate the stochastic optimal control problem 
formulated in the preceding section as one with a fixed planning horizon and to derive a dynamic 
programming principle and the corresponding HJB equation. 

3.1. Dynamic programming principle. Let us denote by A{t,x) the set of admissible de- 
cision strategies i' = {c,p,9) for the dynamics of the wealth process with boundary condition 
X{t) = x. For any u G A{t,x) we define the functional 

U{c{s),s)ds + BiZiT),T)I{^^T} + W{XiT))I{,:,T} r>t,J^t 



J{t,x;u) = Et^x 
and introduce the associated value function 



V{t,x)= sup J{t,x;v) . 

ueAit,x) 

Given some initial condition (t, x) G [0, T] x M, we have that i/* G A{t, x) is an optimal control 
iiV{t,x) = J{t,x;iy*). 

The following lemma is the key tool to restating the control problem above as an equivalent 
one with a fixed planning horizon. See |10j for a proof. 

Lemma 3.1. Suppose that the utility function U is either nonnegative or nonpositive. If the 
random variable t is independent of the filtration ¥, then 



J{t,x;u) = Et^x 



T 



F{s, t)U{c{s), s) + fis, t)B{Z{s), s) ds + F{T, t)W{X{T)) 



where F(s,t) is the conditional probability for the wage earner's death to occur at time s condi- 
tional upon the wage earner being alive at time t < s and f{s,t) is the corresponding conditional 
probability density function. 

Using the previous lemma, one can state the following dynamic programming principle, ob- 
taining a recursive relationship for the maximum expected utility as a function of the wage 
earner's age and his wealth at that time. See [lOj for a proof. 
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Lemma 3.2 (Dynamic programming principle). For < t < s < T, the maximum expected 
utility V{t, x) satisfies the recursive relation 



V{t,x) 



sup E 



A(n) dn V{s,X{s)) 



exp 



+ j F{u,t)U{c{u),u) + f{u,t)B{Z{u),u) du 



The transformation to a fixed planning horizon can then be given the following interpretation: 
a wage earner facing unpredictable death acts as if he will live until time T, but with a subjective 
rate of time preferences equal to his "force of mortality" for the consumption of his family and 
his terminal wealth. 



3.2. Hamilton-Jacobi-Bellman equation. The dynamic programming principle enables us 
to state the HJB equation, a second-order partial differential equation whose "solution" is the 
value function of the optimal control problem under consideration here. The techniques used 
in the derivation of the HJB equation and the proof of the next theorem follow closely those in 

El [ini [IB- 
Theorem 3.3. Suppose that the maximum expected utility V is of class C^. Then V satisfies 
the Hamilton-Jacobi-Bellman equation 

Vt{t,x) - X{t)V{t,x) + sup 'H{t,x;i^) = 

u€A{t,x) , (7) 

V{T,x) = W{x) 
where the Hamiltonian function T-L is given by 

n{t,x;u) = (^i{t)-c-p+ (^r{t) + ^^0nMt)-r{t))jx^V^{t,x) 

2 M / N \2 
+ yY1 [j2Gn<Tnm{t)j V^^it, x) + Xit) B (^X + , tj 



+ U{c,t) . 



Moreover, an admissible strategy u* = {c*,p*,9*) whose corresponding wealth is X* is optimal 
if and only if for a.e. s G [t^T] and P-a.s. we have 

Vt{s,X*{s)) - X{s)V{s,X*{s))+'H{s,X*{sy,iJ*) = . (8) 

Proof. We divide the proof in two parts: we start by establishing the HJB equation ([7]) and then 
we will prove that the equality ([S]) holds. 

Recall that the wealth process X(t), t € [0, min{T, r}], satisfies the stochastic differential 
equation ([S]). Using Ito's lemma, we obtain that 

l-t+h M „t+/i 

Vit + h,X{t + h)) = V{t,X{t))+ a{u,X{u)) du+y^ b^iu,X{u))dWmiu) , {9) 

m=i 
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where the integrand functions a and bm, m = 1, . . . , M, are given by 

a{t,X{t)) = Vt{t,X{t))+ (^i{t)-c-p + X{t) (^r{t) + Y^^6n{fin{t) - r{t))yj VS,X{t)) 

M / N \ 2 

+ -X\t)V^^it,X{t)) 



(10) 



m=l \n=l 



N 



h^{t,X{t)) = X{t)V,{t,X{t))Y,&n(Tnm{t) , m = l,...,M. 



n=l 



Using the dynamic programming principle of Lemma l3.2l and setting s = t + h we get the identity 



V{t, x) = sup E 



exp 



X{u) du V{t + h,X{t + h)) 



t+h 



+ J F{u, t)U{c{u),u) + f{u, t)B{Z{u),u) du 

Noting that for small enough values of h the following inequalities hold 

rt+h 



(11) 



exp ^ \{v) d^^\- X{t)h ± 0{h^) 



[12) 



and combining the equality (jlip with the inequalities ()12p above, we obtain 



€ sup E 

t+h 



(1 - X{t)h ± 0{h^))V{t + h, X{t + h)) - Vit, x) 



G sup E 

u<^A(t,x) 



+ J F{u,t)U{c{u),u) + f{u,t)B{Z{u),u) du 
Substituting ([9]) in the last equation we obtain 

[1 - X{t)h ± 0{h^)) {v{t, X{t)) + a{u, X{u)) du 

M r-t+h 

+ (1 - X{t)h ± 0{h^)) / hraiu, X{u)) dWm{u) 

m=l 

rt+h 

-V{t, x) + J (1 - F{u, t))U{c{u),u) + A(t)(l - Fiu, t))B{Z{u),u) du 
Dividing the previous equation by h and letting h go to zero we obtain the equality 



J't 



= sup 

u&A(t,x) 



Vt{t, x) - X{t)V{t, x) + (^{t) -C-P+ ^r(t) + J2 0nMt) - r{t))j x^ V^{t, x) 

2 M / N \ 2 

+ yY.[Y1 ^nCTnmit) J V^^it, x) + X{t)B {Z{t),t) + U (c, t) 
rn=l \?i=l / 
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Letting Z{t) = x + the dynamic programming equation becomes 



sup 



Vt{t,x) - X{t)V{t,x)+ (^i{t)-c-p + (^r{t) + Y^^9n{fin{t) - r{t))j V^{t, 

2 M / N \2 
+ yY1 ( V^^it, x) + X{t) B (^X + , tj + U {c, t) 

Finally, noting that Vt{t,x) — X{t)V{t,x) does not depend on z^, we obtain the HJB equation of 
([7|), thus concluding the proof of the first part of the theorem. 

Regarding the second part of the theorem, we start by letting G A{t, x) and apply Ito's 
lemma to 



exp 



X{v) dv V{s,X{s)) . 



We obtain that 
V{t,x) 



exp 



exp 



X{v) dv\w{X{T)) 



X{v) dv {a{u,X{u)) - X{u)V{u,X{u))) du 



■E 

m=l 



exp 



X{v) dv bm{u,X{u)) dWm{u) 



where a{t,x) and bmit,x), m = 1, . . . ,M are as given in (fTO]l . From the previous equality, we 
get 

V{t,x) 



t,x 



exp 

T 

exp 



X{v) dvj W{X{T)) 
X{v) dv] {a{u,X(u)) - X{u)V{u,Xiu))) du 



(13) 



From the definition of the hazard function in ([2|), we obtain that the conditional probability 
-F(s, t) for the wage earner's death to occur at time s conditional upon the wage earner being 
alive at time t < s is given by 



F(s, t) = exp 



X{v) dv 



(14) 



Similarly, we obtain that the conditional probability density function f{s,t) for the death to 
occur at time s conditional upon the wage earner being alive at time t < s is given by 



f{s,t) = X{s)exp {- j'^ 



X{v) dv 



Using (fn|) and we rewrite (fTHj) as 



F{T,t)W{X{T))- / F{u,t){a{u,X{u)) - X{u)V{u,X{u))) du 



V{t,x) = Et,^ 



Using Lemma |3.H we rearrange (|16p to obtain 



(15) 



(16) 



V{t,x) = J{t,x;u) 
t-T 



't,x 



F{u, t) {Vt{u, X{u)) - X{u)V{u, X{u)) + n{u, X{u);u)) du 



(17) 
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Consider now an optimal admissible strategy z/* = (c*,p*,9*) whose corresponding wealth is 
X* . From equation ()17p we have that 

V{t,x) = J{t,x;v*) 



(18) 



F{u,t){Vt{u,X*{u))-\{u)V{u,X*{u))^n{u,X*{u\,v*)) dn 

and from the HJB equation d?]), we have 

Vt(u, X*{u)) - \{u)V{u, X*{u)) + U{u, X*(u); v*) < 0. (19) 

Combining (jlSp and ()19p , we obtain that u* is optimal if and only if the value function V satisfies 
([8]), which concludes the proof. □ 

The second part of the theorem above provides a clear approach for the computation of 
optimal insurance, portfolio and consumption strategies. In particular, we obtain the existence 
of such optimal strategies under rather weak conditions on the utility functions. 

Corollary 3.4. Suppose that the maximum expected utility V is of class and that the utility 
functions U and B are strictly concave with respect to their first variable. Then the Hamiltonian 
function % has a regular interior maximum u* = {c*,p*,9*) G A{t,x). 

Proof. Using the second part of theorem (|3.3p . an optimal admissible strategy v* = {c*,p*,9*) 
with wealth process X* must satisfy ([S]). Therefore, i^* must be such that H attains its maximum 
value. We start by remarking that the condition to obtain the maximum for T-L decouples into 
three independent conditions, as seen in the following: 

supn{t,x;i/) = {r{t)x + i{t))Vx{t,x) + sup I^U{c,t) - cVr,{t,x)^ 

+sup |A(t)i3(^x + ^,t^ -pl4(t,x)| (20) 



p 



2 M / N \ 2 N 



+SUp j Y X] ^n(Tnm{t) V^ccit, x) + ^ O^Mt) " r{t))xVx{t, x) \ . 

^ ^ m=l \n=l J n=l ' 

Therefore, it is enough to study the variation of Ti with respect to each one of the variables c, p 
and 6 independently. Thus, computing the first-order conditions for a regular interior maximum 
of 7i with respect to c, p and 9 we obtain, respectively, the following three conditions 

-Vx{t,x) + Ucic\t) = 
-Vx{t,x) + ^Bz(x + ^,t) = (21) 



r]{t) "V Vit)' 
xVx{t, x)a + x^Vxx{t, x)aa^O* 

where the subscripts in U and B denote differentiation with respect to each function's first 
variable, a denotes the risk premium ([T|), and O^n denotes the origin of M.^ . Computing the 
second derivative with respect to each variable (or the Hessian matrix in the case of 9), we 
obtain 

7icc{t,x;u*) = Ucc{c*,t) 

^pp(*'^;^*) = ^fzz{x + ^yt^ (22) 

'Heeit,x;i'*) = x'^Vxxit,x)aa'^ . 

Note that T-iccit^x^v*) is negative since U is strictly concave on its first variable and that 
Tippit, x; I/*) is negative since A(t) is positive for every < t <T and B is strictly concave on its 
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first variable. To see that T-Lo0{t, x; u*) is negative definite, recall that aa'^ is assumed to be non- 
singular and, thus, positive definite. Moreover, note that Vxxit,x) must be negative: if Vxx(t,x) 
is positive, then T-L would not be bounded above and as a consequence of the HJB equation 
either Vt{t,x) or V{t,x) would have to be infinity, contradicting the smoothness assumption on 
V. Therefore, Tigg is negative definite and H has a regular interior maximum. □ 

4. The family of discounted CRRA utilities 

In this section we describe the special case where the wage earner has the same discounted 
CRRA utility functions for the consumption of his family, the size of his legacy, and the size of 
his terminal wealth, i.e., from now on we assume that the utility functions are given by 

U{c, t) = e-P'— , BiZ, t) = e-P'— , W{X) = e""^ , (23) 

7 7 7 

where the risk aversion parameter 7 is such that 7 < 1, 7 7^ 0, and the discount rate p is positive. 

4.1. The optimal strategies. Using the optimality criteria provided in theorem l3.3l we obtain 
the following optimal strategies for discounted CRRA utility functions. 

Proposition 4.1. Let ^ denote the non-singular square matrix given by {aa^)^^. The optimal 
strategies in the case of discounted constant relative risk aversion utility functions are given by 

c*{t,x) = ^{x + hit)) 

p*{t,x) = r?(t)(p(t)-l)x + D(t)6(t)) 

e*{t,x) = ^ {x + b{t))Ca{t) , 
x(l -7) 



where 



rT 

b{t) = I z(s)exp(— / r(v) + rj{v) dv ] ds 



m 

e{t) 



1 (\{t) 



1/(1-7) 



exp H{v) d^ + exp {- H{v) dv^ K{s) ds 



m = a^mait)-^\\a^Ca{t)f . 

Proof. Assume that the utility functions U, B and W are as given in (|23|) . Using the first order 
conditions in (pT|) . we obtain that the optimal strategies depending on the value function V are 
given by 

c*{t,x) = {eP'Vx{t,x))-^'^^-^^ 

^ ,,)((*)^)-""-"-.) ,24) 
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We are now going to find an explicit solution for the HJB equation ([7]). We substitute c, p 
and 9 in the HJB equation by the optimal strategies in (|24p and combine similar terms to arrive 
at the following partial differential equation 

Vt{t, x) - X{t)V{t, x) + ((r(t) + 7]{t))x + i{t)) V^{t, x) 

-m^^^j^ + i^e-^*/(i-^)i^(t)(y.(t,x))-^/(^-^) = , (25) 

where and K[t) are as given in the statement of this proposition and the terminal condition 
is given by 

V{T,x) = W{x). (26) 

We consider an ansatz of the form 

V{t,x) = '^{x + b(t)r , (27) 
7 

and substitute it in (|25p so that a{t) and b[t) are determined by the differential equation 

1 da(t) a{t) dh{t) _ ^.^^a(^) Mt) + i^{t))x + i{t)]a{t) 
7 dt x + h{t) dt ^ ' x + h{t) 

+S(t)^ + i^e-''*/(i-^)i^(t)(a(t))-^/(i-^) = . 
1-7 7 

Note now that the previous differential equation and the terminal condition ()26p decouples into 
two independent boundary value problems for a{t) and h{t) which are given, respectively, by 

+ Ut) + ,(t) _m + m\ ait) + i^e-*/(^--)K(t)(a(t))-/(^-) = 
a{T) = e-P^ , (28) 

and 

- {r{t) + rj{t))h{t) + i{t) = Q 
h{T) = . (29) 
To find a solution for the boundary value problem ()28p . we write a{t) in the form 

a{t) = e-P\e{t)Y~^ , 
obtaining a new boundary value problem for the function e{t) of the form 

^-H{t)e{t)+K{t) = Q 

e{T) = 1 , (30) 

where K{t) and H[t) are as given in the statement of this proposition. Since equation ()30p is 
a linear, non- autonomous, first order ordinary differential equation, it clearly has an explicit 
solution of the form 



e{t) = exp H{v) di;^ + ^ exp H{ 



(v) dv K{s) ds 



Therefore, we obtain that the solution of (|28p is given by 

a{t) = e-P^ (^exp H{v) d^ + ^ exp H{v) dv^ K{s) ds^ . (31) 
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To find a solution for the boundary value problem ()29p . we just note that this is again a linear, 
non-autonomous, first order differential equation and its solution is given by 

b{t) = i{s) exp ^- r{v) + ■q{v) dv^ ds (32) 

as required. 

Combining (|24p with (|27p . (|3ip and (j32p . we obtain that the optimal strategies in the case of 
CRRA utilities are then given by 

c*{t,x) = -l-^{x + b{t)) 

p*{t,x) = r,{t)i{Dit) - l)x + Dit)b{t)) 

e*{t,x) = ^±^Cait), 

where D{t) is as given in the statement of this proposition, which concludes the proof. □ 

Note that the quantities b{t) and x + b{t) are of essential relevance for the definition of the 
optimal strategies in proposition 14.11 The quantity b(t), that we will refer to as human capital 
following the nomenclature introduced in [B], should be seen as representing the fair value at 
time t of the wage earner's future income from time t to time T, while the quantity x + b{t) 
should be thought of as the full wealth (present wealth plus future income) of the wage earner 
at time t. It is then natural that these two quantities play a central role in the choice of optimal 
strategies, since they determine the present and future wealth available for the wage earner and 
his family. 

Remark 4.2. Noting that r{t), r]{t) and i{t) are positive functions and considering the boundary 
value problem (|32p . if r{t) + rj{t) is small enough, we can deduce that human capital function 
b{t) has the following properties: 

a) it is a positive function for all < t < T; 

b) it is concave. 

Moreover, we have that b{t) is either: 

i) a decreasing function for all t E [0, T]; or 

ii) a unimodal map oft, i.e. there exists some t* € (0,r) such that b{t) is increasing for 
all < t < t* , decreasing for all t* < t < T; Furthermore, we have that the graph of b(t) 
intersects the graph of the function i(t)/(r{t) + ri{t)) at t = t* . 

From the explicit knowledge of the optimal strategies, several economically relevant conclu- 
sions can be obtained. See Fig. [T] for a graphical representation of the optimal life-insurance 
purchase as a function of age and "full wealth" x + b{t) of the wage earner. We start by prov- 
ing an auxiliary lemma before moving on to the statement and proof of some of the optimal 
strategies properties. 

Lemma 4.3. Suppose that for all t E [0, min{T, r}] the following two conditions are satisfied: 

a) \{t) < r](t); 

b) H{t) < 1. 

Then, the inequality D{t) < 1 holds for every t E [0, min{T, r}]. 
Proof. Recall that D[t) is given by 

1 [\{t)V'^'-^^ 




Figure 1. The optimal life-insurance purchase for a wage earner that starts 
working at age 25 and retires 40 years later. The parameters of the model were 
taken as iV = M = 2, i{t) = 50000 exp(0.03t), r = 0.04, p = 0.03, 7 = -3, 
X(t) = 0.001+exp(-9.5 + 0.1i), r]{t) = 1.05A(t), m = 0.07, ^2 = 0.11, au = 0.19, 
CT12 = 0.15, (T21 = 0.17 and 1T22 = 0.21. 



Using condition b) and noting that K{t) is positive for all t € [0, min{T, r}], we have that 
e{t) = exp H{v) dv^ + ^ exp H{v) dv^ K{s) ds 

> exp ^- ^ Idv^ + exp ^ 1 K{s) ds > 1 . 

Putting together the previous inequality and condition a) , we obtain the required inequality. □ 

The next result provides a qualitative characterization of the optimal life insurance purchase 
strategy. 

Corollary 4.4. Assume that the conditions of lemma \4-3\ are satisfied. Then, the optimal 
insurance purchase strategy p*{t,x) has the following properties: 

a) it is a decreasing function of the wealth x; 

b) it is an increasing function of the wage earner's human capital b{t); 

c) it is negative for suitable pairs of wealth x and "age" t; 

d) if the wage earner's wealth x is small enough and r]{t) is non- decreasing, the function 
t I— 7- p*{t,x — b{t)) has the same monotonicity as the human capital function b{t). 

Proof. Recall from proposition 14.11 that the optimal insurance purchase strategy p*{t, x) is given 
by 

p*{t, x) = r,{t) {{D{t) -l)x + D{t)b{t)) (33) 

for all t G [0,min{T,r}]. 
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Items a) and b) follow from lemma since D(t) is a positive function such that D{t) < 1 
for all t G [0,min{r,r}]. 

For the proof of item c), note that p*{t,x) is negative for all {t,x) E [0,T] x such that 

^ ' Mt) > , 



e(t)r/(t)i/(i-7) _ A(t)V(i-7) 
and positive otherwise. 

Item d) follows from ()33p and the fact that ri{t) is a non-decreasing. □ 

Some comments regarding the assumptions in lemma 14.31 (and corollary 14. 4p seem necessary. 
Starting with condition (a), the life insurance company must establish the premium-insurance 
ri{t) in such a way that A(t) < ry(t) in order to make a profit (the insurance policy being fair 
whenever X{t) < r]{t)). Regarding condition (b), we note that the quantities r, p, rj and A are 
usually very small in the real world and, moreover, the relative risk aversion of the wage earner 
is negative in general. This is consistent with the assumption that H{t) is bounded above by 
some positive constant. 

Apart from studying how optimal life insurance purchase varies with age and wealth, it is also 
relevant to understand how the remaining parameters which define the financial and insurance 
markets influence life insurance purchase. 

Corollary 4.5. With all other parameters, including t andx constant, the optimal life insurance 
purchase rate p*{t,x) is an increasing function of the discount rate p. 

Proof. To study the influence of the discount rate p on the optimal optimal life insurance pur- 
chase, we consider two different values of p and compare the corresponding values of the optimal 
life insurance purchase. We distinguish the functions associated with each of the two parameter 
values by their subscript. 

Assume that pi and p2 are such that pi < p2. Recall the definitions of p*(t, x), H{t), e{t) and 
D(t) given in proposition 14.11 Then, it is clear that the inequalities 

Hi{t) < H2{t) 
eiit) > e2{t) 
Di{t) < D2{t) 

hold for all t G [0, min{T, r}], where the subscripts correspond to pi and p2 in an obvious 
manner. Rewriting x) as 

p*it,x) = 7]it){Dit){x + bit))-x) , 
it follows from the preceding inequalities that 

Pl{t,x) < p*2{t,x) 

for all t € [0, min{T, r}], concluding the proof of the statement. □ 

Remark 4.6. The variation of the optimal life insurance purchase rate p*{t,x) with respect 
to the interest rate r{t), the risk aversion parameter 7, the hazard rate \{t) and the insurance 
premium-payout ratio ri{t) is non-trivial. However, by studying the function p*{t,x) given in 
proposition \4.1\ we can make the following observations: 

i) p*{t,x) is a decreasing function of the interest rate r[t), except for large values of x and 
t close enough to T; 

ii) p*{t,x) is a decreasing function of the risk aversion parameter 7, except for values oft 
close enough to T; 
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iii) p*{t,x) is an increasing function of the hazard rate X{t) and the insurance premium- 
payout ratio rj(t) for small enough values of wealth x and a decreasing function for large 
values of X. 

The extra risky securities in our model introduce novel features to the wage earner's portfolio 
management, as is exemplified in the following result. 

Corollary 4.7. Let ^ denote the non-singular square matrix given by {^crcr^^ ^ and let (^Qi(t))j^ 
denote the n-th component of the vector ^a{t) . The optimal portfolio process 9* {t,x) = (0*,...,^^) 
is such that for every n E {1, A^}." 

a) 0* has the same sign as {^a{t))n; 

b) 0* is a decreasing function of the total wealth x if > for all t G [0,min{r, r}] 
and an increasing function of x if {(,a{t))n < for all t E [0, minjT, r}]; 

c) 0* is an increasing function of the wage earner's human capital b{t) if {(,a{t))n > 
for all t G [0, min{r, r}] and a decreasing function of b{t) if (^a(t))„ < for all t G 
[0,min{r,T}]. 

Furthermore, for every n,m {!,..., A'^} the following equalities hold 

lime:it,x) = +00 1^ = 14% 

x^o+ ^ x^o+9:^{t,x) {^a{t))m 

lime:{t,x) = IimOt,x) = M^. 

a;^oo 1-7 I-7 

Proof. Items a), b) and c) on the first part of the corollary follow from the form of 0*, n € 
{1, A^}, given in the statement of proposition 14.11 and positivity of b{t). 

The limiting behaviours on the second part of the corollary also follow from the form of 0*, 
n£{l,...,N}. □ 

Remark 4.8. Corollary \4. 7| is a mutual fund result: the relative proportions among the risky 
securities are independent of all parameters except for the interest rate and the risky assets 
appreciation rates and volatilities since, for any n,m {1,..., A^} we have 

0*n{t.x) ^ {ia{t))n 

e*^it,x) i^ait)u ■ 

The quantities {(,a{t))n, n € {1,...,A^} can be though as "weighted risk premiums" for the 
risky assets Si, ... , Sn, where the weights are provided by (quadratic) functions on the coeffi- 
cients of the matrix of risky- assets volatilities a. 

Under the assumption that the "weighted risk premiums" {^a{t))n are positive for all n € 
{1, A^} and t € [0, minjT, r}], we obtain the following interesting consequence of the previous 
corollary. 

Corollary 4.9. Let ^ denote the non-singular square matrix given by {aa'^)~^ and let {S,a{t))n 
denote the n-th component of the vector ^a{t). Assume that for every n G {1,..., A^} we have 
{^a{t))n > for all t G [0, min{T, r}] . Then the optimal strategy for wage earners with small 
enough wealth x is to short the risk-free security and hold higher amounts of risky assets. 

Proof. The corollary follows from corollary 14.71 and the fact that the wage earner has no budget 
limitations, thus allowing him to get into short positions on the risk free asset 5*0 of arbitrary 
size. □ 

We conclude this session with a result concerning some qualitative properties of the optimal 
consumption strategy. Its proof follows trivially from the form of the optimal consumption 
c*{t,x) given in proposition 14.11 

Corollary 4.10. The optimal consumption rate c*{t) is an increasing function of both the wealth 
X and the human capital b{t). 



OPTIMAL LIFE INSURANCE PURCHASE 



17 



4.2. The interaction between life insurance purchase and portfoUo management. In 

this section we compare the optimal hfe-insurance strategies for a wage earner who faces the 
following two situations: 

a) in the first case, we assume that the wage earner has access to an insurance market 
as described above and that his goal is to maximize the combined utility of his family 
consumption for all t < minjT, r}, his wealth at retirement date T if he lives that long, 
and the value of his estate in the event of premature death. The optimal strategies for 
the wage earner in this setting are given in Proposition 14.11 

b) in the second case, we assume that the wage earner is without the opportunity of buying 
life insurance. His goal is to maximize the combined utility of his family consumption for 
all t < min{r, r} and his wealth at retirement date T if he lives that long. Similarly to 
what we have done previously, we translate this situation to the language of stochastic 
optimal control and derive explicit solutions in the case of discounted CRRA utilities. 

We concentrate on the case b) described above for the moment. Similarly to what was done in 
case a), this problem can be formulated by means of optimal control theory. The wage earner's 
goal is then to maximize a new cost functional subject to: 

• the dynamics of the state variable, i.e., the dynamics of a wealth process X^{t) given by 



where t G [0, minjT, r}] and x is the wage earner's initial wealth. 

• constraints on the remaining control variables, i.e., the consumption process c^{t) and 
the reduced portfolio process 9'^{t) = (6'?(t), • • • , 6i^(t)) G M^; and 

• boundary conditions on the state variables. 

Let us denote by A^{x) the set of all admissible decision strategies, i.e. all admissible choices 
for the control variables = {c^,9^) € M^^^. The dependence of A^{x) on x denotes the 
restriction imposed on the wealth process by the boundary condition X^(0) = x. 

The wage earner's problem can then be stated as follows: find a strategy = {c^, 9^) G A^{x) 
which maximizes the expected utility 



where U{c^, •) is again the utility function describing the wage earner's family preferences re- 
garding consumption in the time interval [0,min{T, r}] and B{X^,t) is the utility function for 
the terminal wealth at time t = T At. As before, we restrict ourselves to the special case where 
the wage earner has the same discounted CRRA utility functions for the consumption of his 
family and his terminal wealth given in (|23p . The optimal strategies are given in the next result 

Proposition 4.11. Let ^ denote the non-singular square matrix given by {aa^)~^ . The optimal 
strategies for problem (|34p in the case where U{c^,-) and B{X^ , ■) are the discounted constant 
relative risk aversion utility functions in ([25]) are given by 





■TAt 



(34) 



c^*{t,x) 




{x + b'{t)) 



9^\t,x) 



1 



{x + b\t))^a{t) 



x(l -7) 
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where 

rT 



b°{t) = i{s)exp r{v) dv^ ds 

e°(t) = exp(^-^ H^{v)dv^+j^ exp (v) dv^ ds 



^ ^ 1-7 (1-7)2 1-7 



and is as given in the statement of proposition 

We skip the proof of the previous proposition, since it is of the same nature as the proofs of 
theorem 13.31 and proposition 14.11 

We should point that it would be preferable from an economic standpoint to maximize also the 
utility of wealth at the time of premature death. In fact, it is easy to check that this corresponds 
to the addition of a constraint of the form p{t, x) = to our original problem. However, the 
corresponding HJB equation would contain an extra term of the form X(t) x'^/j. The presence 
of this additional term in the HJB equation makes it impossible for us to obtain closed form 
solutions. On the other hand, a strategy optimizing the final wealth at retirement should be 
close to an optimal strategy which maximizes also the wealth for the case of an eventual death 
of the wage earner before retirement time. We plan to address such a comparison in future 
research. 

If we make the wage earner income i{t) equal to zero, then we get a solution which is close 
to Merton's classical solution, but it still depends on the hazard function X{t), i.e., even in the 
absence of a life insurance policy the uncertainty regarding wage earner's lifetime still plays a 
role on the determination of the optimal consumption and investment strategies. 

Propositions 14. ll and 14. iT] provide us with optimal portfolio processes for the two settings a) 
and b) described above. In the next theorem we show how these optimal portfolio processes 
compare. 

Theorem 4.12. Let ^ denote the non-singular square matrix given by {aa'^)~^ and (^a(t))„ the 
n-th component of the vector ^a{t). For each n G {1,...,A^}, we have that 6^*{t,x) > 0^{t,x) if 
and only if {^a{t))n > 0. 

Proof. Recall the definitions of 9^* and 9* from propositions 14.111 and 14.11 respectively. Note 
that 

9r{t,x)-9Ut,x) = ^ (^ + 6°(t))(g«(t))„- ^ (x + b(t))(ga(t))„ 

x[l — 7j x{l — 7j 

= |r5^(^°w-^w)- (35) 

Using the definitions of b^{t) and b{t) given in the statements of propositions 14.111 and 14. H 
respectively, we obtain that 

6°(t) - b{t) = [ i(s)exp (- [ r{v) dv^ (l - exp (- [ r]{v) dv]] ds . (36) 



Jt \ Jt / \ 

Since rj{t) is a positive function, we get that 

j r]{v) dv > 

and therefore, the inequality 

1 - exp ( - / r]{v) dt; ) > (37) 



OPTIMAL LIFE INSURANCE PURCHASE 



19 



holds for all < t < s < r. Therefore, combining ()35p . ()36p and ()37p and recalling that 
7 < 1, we obtain that the sign of 9^ it,x) — 6^(t,x) is the same as the sign of {£,a{t))n for all 
t € [0,min{r,T}]. □ 

The economic implications of the theorem above are made clear in the following result. 

Corollary 4.13. Let ^ denote the non-singular square matrix given by {aa^)"^ and (Ca(^))n 
the n-th component of the vector ^a{t). Assume that for every n G {1,...,A^} we have that 
{^a{t))n > 0. Then, the optimal portfolio of a wage earner with the possibility of buying a life 
insurance policy is more conservative than the optimal portfolio of the same wage earner if he 
does not have the opportunity to buy life insurance. 

5. Conclusions 

We have introduced a model for optimal insurance purchase, consumption and investment 
for a wage earner with an uncertain lifetime with an underlying financial market consisting 
of one risk-free security and a fixed number of risky securities with diffusive terms driven by 
multidimensional Brownian motion. 

When we restrict ourselves to the case where the wage earner has the same discounted CRRA 
utility functions for the consumption of his family, the size of his legacy and his terminal wealth, 
we obtain explicit optimal strategies and describe new properties of these optimal strategies. 
Namely, we obtain economically relevant conclusions such as: (i) a young wage earner with 
smaller wealth has an optimal portfolio with larger values of volatility and higher expected 
returns; and (ii) a wage earner who can buy life insurance policies will choose a more conservative 
portfolio than a similar wage earner who is without the opportunity to buy life insurance. 

It is worth noting that the model described in this paper can be improved by adding further 
ingredients such as, for instance, some form of stochasticity on the income function or the 
hazard rate, some constraints on the possibility of trading of stocks on margin, the existence of 
additional life insurance products on the market or even adding some correlation between the 
wage earner's mortality rate and the underlying financial market. We plan to address at least 
some of these issues in a future publication. 
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